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"j^ ' We investigate a distributed optimal control problem for a phase field 

model of Cahn-Hilliard type. The model describes two-species phase 
segregation on an atomic lattice under the presence of diffusion; it has 
been introduced recently in [3] , on the basis of the theory developed in 
J> ■ |15| . and consists of a system of two highly nonlinear ly coupled PDEs. 

00 i For this reason, standard arguments of optimal control theory do not 

| apply directly, although the control constraints and the cost functional 

£f~) • are of standard type. We show that the problem admits a solution, 

\q ■ and we derive the first-order necessary conditions of optimality. 

o 

1 Introduction 

• i-H , 

K> . 

Let Q C IR denote an open and bounded domain whose smooth boundary 
T has outward unit normal n, let T > be a given final time, and let 
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Q := Q x (0, T), E := T x (0, T) . In this paper, we study distributed 
optimal control problems of the following form: 

(CP) Minimize the cost functional 

J(u,p,p) = - I \p(x,T) - p T (x)\ 2 dx 
z Jn 

+ / / \p(x, t) — pr(x, t)\ 2 dx dt 

2 Jo Jn 

+ 1T f I Hx,t)\ 2 dxdt (1.1) 

2 Jo Jn 

subject to the state system 

(e + 2 p)p t + pp t - Ap = u 
Spt -Ap + f{p) = p 

^1 = ^ = 
dn dn 

p(x, 0) = p (x) , fi(x,0) = p (x) , 
and to the box control constraints 

u e U ad = {u e L°°(Q) ; < u < U a. e. in Q} . (1.6) 

Here, /?i>0, /3 2 >0, e>0, and 5 > are constants; U G L°°(Q) de- 
notes a given bound, and G L 2 (f2) and /iy G L 2 (Q) represent prescribed 
target functions of the tracking- type functional J . Although for large parts 
of the subsequent analysis much more general cost functionals could be ad- 
mitted, we restrict ourselves to the above situation for the sake of a simpler 
exposition. 

The state system (jl.2ft - fll.5p constitutes a phase field model of Cahn-Hilliard 
type that describes phase segregation of two species (atoms and vacancies, 
say) on a lattice in the presence of diffusion; it has been introduced recently 
in [15, 4]. The state variables are the order parameter p, interpreted as a 
volumetric density, and the chemical potential p . For physical reasons, we 
must have < p < 1 and p > almost everywhere in Q . The control 
function u on the right-hand side of ( 11. 2p plays the role of a microenergy 
source (see below). Moreover, the nonlinearity / is a double- well potential 
defined in (0,1), whose derivative /' is singular at the endpoints p = and 
P = 1: e.g., / = /i + f 2 , with f 2 smooth and f x {p) = c(p log(p) + (1 - 
p) log(l — p)) , with c a positive constant. 



a.e. in Q, (1.2) 

a. e. in Q, (1.3) 

a. e. on E, (1.4) 

a.e. in Q, (1-5) 
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System (ll.2l) - (ll.5p is singular, with highly nonlinear and nonstandard cou- 
pling. In particular, nasty nonlinear terms involving time derivatives occur 
in ( II. 2p . and the expression f'(p) in (II. 3p may become singular. For the case 
u = (no control), this system was analyzed in a recent paper pQ; the case 
£ \ was studied in [5]. We also refer to the papers [2] and [3], where the 
corresponding Allen-Cahn model was discussed. 

The mathematical literature on control problems for phase field systems is 
scarce and usually restricted to the so-called Caginalp model of phase tran- 
sitions (see, e.g., [11], [9], [10], [17], and the references given there). More 
general, thermodynamically consistent phase field models were the subject 
of [13]. Control problems for the system (11.2p - (11.5p have never been stud- 
ied before. We remark at this place that it would be a challenging task to 
study boundary control problems for the PDE system (11.21) . (II. 3p in place of 
distributed ones as in this paper; notice, however, that this would require to 
first establish appropriate well-posedness results for non-homogeneous Neu- 
mann boundary conditions or for non- homogeneous boundary conditions of 
third kind. Such results are presently not available. 

The paper is organized as follows: below, we briefly recall the thermodynamic 
background of the state system (ll.2l) - (ll.5p . In Section 2, we establish the 
existence of a solution to the optimal control problem. First-order necessary 
optimality conditions, as usual given in terms of the adjoint system and a 
variational inequality, are derived in Section 3. A large part of this analy- 
sis is devoted to proving that the control-to-state mapping is directionally 
differentiable in appropriate function spaces. 

1.1 Some thermodynamic background 

The state equations (II. 2p . (11.31) result from the balances of microenergies 
and microforces postulated in a model for phase segregation and diffusion of 
atomic species on a lattice introduced in [T5] , a paper we refer the reader to 
for details. That model is a variation of the Cahn-Hilliard system 



when, for the sake of simplicity, the mobility coefficient k > is taken equal 
to one. Customarily, the equations in ( II. 7p are combined so as to get the 
well-known Cahn-Hilliard equation 



which describes diffusive phase separation processes in a two-phase material 
body. 



Pt-K Ap = 0, p = -Ap + f'(p), 
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A generalization of (jl.8p was introduced by Fried and Gurtin in the papers 
[B] and [S]. Here is their line of reasoning: 

(i) to regard the second equation in (11. 7p as a balance of microforces : 

div£ + 7r + 7 = 0, (1.9) 

where the distance microforce per unit volume is split into an internal part 
7r and an external part 7, and where the contact microforce per unit area 
of a surface oriented by its normal n is measured by £ ■ n in terms of the 
microstress vector 

(ii) to interpret the first equation in (ll.7p as a balance law for the order 
parameter: 

p t = -divh + a, (1.10) 
where the pair (h,o~) is the inflow of p; 

(iii) to restrict the admissible constitutive choices for n,£,h, and the free 
energy density if), to those consistent in the sense of Coleman and Noll [I], 
with an ad hoc version of the Second Law of Thermodynamics - namely, a 
postulated "dissipation inequality that accommodates diffusion" - given in 
the form 

4>t + (n - ti) Pt - £ ■ Vpt + h- Vp< (1.11) 

(cf., in particular, Eq. (3.6) of [8]). Within this framework, an admissible 
set of constitutive prescriptions turns out to be: 

ip = $(p,Vp), 7r(p, Vp,p) = fi - d p $(p, Vp), £(p, Vp) = 9 Vp ^(p,Vp), 

(1.12) 

together with 

h = -MVfi, where M = M(p, Vp, p, Vp). (1.13) 

Moreover, it follows that the tensor-valued mobility mapping M must obey 
the inequality 

Vp • M(p, Vp, p, Vp) Vp > 0. 
It follows from (Ojl . fOO]) . (TTTT2D . and ffTT3|) 1 that 

Pt = div (MV (^(p, Vp) - div(<9 V/ ^(p, Vp)) - 7)) + <r ; 
the Cahn-Hilliard equation (11. 8p results for the special choice 

^(p,Vp) = /(p) + ^|Vp| 2 , M = «l, (1.14) 
1 



provided that the external distance microforce 7 and the order parameter 
source term a are taken identically zero. 

In contrast to the theory developed by Fried and Gurtin, the approach taken 
in [15] was the following: while step (i) was retained, the order parameter bal- 
ance fll.lOp and the dissipation inequality (II. lip were replaced, respectively, 
by the microenergy balance 

e t = e + w, e := — div h + a, w := —n p t + £ • V p t , (1-15) 

and the microentropy imbalance 

f] t > — div h + er, h := ph, a := pa. (1-16) 

The salient new feature of this approach to phase segregation modeling is that 
the microentropy inflow (h, a) is deemed proportional to the microenergy 
inflow (h,a) through the chemical potential p; consistently, the free energy 
is defined to be 

i/j:=e-p-\ (1.17) 

where the chemical potential plays the same role as the coldness in the deduc- 
tion of the heat equation. Just as the absolute temperature is a macroscopic 
measure of microscopic agitation, its inverse - the coldness - measures micro- 
scopic quiet. Likewise, as argued in [15], the chemical potential can be seen 
as a macroscopic measure of microscopic organization; and, just as is always 
done for coldness, one can provisionally assume that p is positive almost 
everywhere in Q . This assumption, which is important to proving that the 
resulting system of field equations does have solutions, must be justified a 
posteriori. The requirement that p be positive is also the reason why we 
cannot admit negative controls u in the control problem (ll.ip - (11.6p . 

Combining (II . 15p -( TT . 1 7[) . and assuming that p > 0, one finds that 

^t< -v dtip- 1 ) + p^h ■ Vp-Txpt + Z- V Pt ; (1.18) 

this reduced dissipation inequality replaces (jl.lip in filtering out a la Coleman- 
Noll the inadmissible constitutive choices. 

On taking all of the constitutive mappings delivering ir,£, 77, and h, to 
depend in principle on the list of variables p, Vp, p, V/i, and on choosing 

V = f (p, Vp,p) = -p P + f(p) + ||Vp| 2 , (1.19) 
one sees that compatibility with ( I1.18P implies that 

vr(p, Vp, p) = -d p t/j(p, Vp, p) = p - f'(p), 

?(P> V P> A 1 ) = ^v P ^(p, Vp, p) = Vp, 

V(P, Vp, p) = p 2 d^(p, Vp, p) = -p 2 p, (1-20) 
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together with 



h(p, Vp, p, Vp) = H(p, Vp, p, Vp) Vp, Vp • H(p, Vp, p, Vp) Vp > 0. 

If we now choose for H the simplest expression H = nl , implying a con- 
stant and isotropic mobility, and if we once again assume that the external 
distance microforce 7 and the source a are null, then we can infer from 
(ll.20p and (I1.17P that the microforce balance (11. 9p and the energy balance 
(I1.15P become, respectively, 

div(Vp) + p-/'(p)=0, (1.21) 
2pp t + pp t - nAp = 0. (1.22) 

This is a nonlinear system for the unknowns p and p , to be compared with 
system (II. 7ft : while equations (ll.2ip and (II. 7p 9 coincide, equation (ll.22p is 
considerably more difficult to handle than (11.71) 1 . Indeed, the latter is linear 
while the former is not; moreover, the time derivatives of p and p are both 
present in (ll.22p . and there are nonconstant factors in front of both p, t and 
p t that should remain positive during the entire evolution. Note that, for 
nonzero microenergy source a, Eq. (11.221) becomes: 

2 p pi + p p t — n Ap = —a. (1-23) 

In this sense, the control variable u in (ll.2p is nothing but —a . 
So far, it has not been possible to tackle the system (11.211) . (11.221) (nor (11.211) . 
(I1.23P ) mathematically. Not so for system (11.21) . (11.31) . a regularized version 
of (11.211) . (I1.23P (with u = —a ) obtained by introducing the extra terms e d t p 
in (11.231) and 5 dtp in (ll.2ip . with small positive coefficients e and 5 (our 
motivations for including such terms have been proposed and emphasized in 

my 

2 Problem statement and existence 

Consider the optimal control problem fl 1 . 2 j) — ( TTT6]) . For convenience, we intro- 
duce the abbreviated notation H = L 2 (Q), V = H\tt) , W = {w E H 2 (tt) ; 
dw/dn = on T} . We endow these spaces with their standard norms, for 
which we use self-explaining notation like || • ||y ; for simplicity, we also write 
|| ■ \\h for the norm in the space H x H x H . Recall that the embed- 
dings W C V C H are compact. Moreover, since V is dense in H , we 
can identify H with a subspace of V* in the usual way, i.e., by setting 
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(u,v)v*,v = ( u j v )h for all u E H and v E V, where (• , •)v*,v denotes 
the duality pairing between V* and V . Then also the embedding H C V* 
is compact, and since iV < 3, we have the continuous Sobolev embeddings 
W C C(H) and V C L 6 (fi) . 

We make the following assumptions on the data: 

(Al) f = fi + f 2 , where f\ E C 2 (0, 1) is convex, f 2 E C 2 [0, 1] , and 

Kmf[(r) = -oo, lim /{(r) = +oo. (2.1) 

(A2) po e W, f'(po) e H, p E V f) L°°(Q) , and 

inf {po(x); x G f2} > 0, sup{po(x); x G fl} < 1, po > a. e. in Q. 

(2.2) 

Notice that (A2) implies that po G C(f2), and that the convexity of fx 
implies that f(po) E H . 

An argumentation that parallels (and thus needs no repetition) the lines of 
the proofs of Theorem 2.2 and Theorem 2.3 of [I] (where we had u = 0) 
shows that the following well-posedness result holds for the state system 
(TOHT3]): 

Theorem 2.1 Suppose that the hypotheses (Al) and (A2) are satisfied. 
Then we have: 

(i) For every u E U a( i the state system (jl.2p - (|1.5|) has a unique solution 
(p, p) such that 

p E W llOC (0, T; if) n ^(0, T; V) PI L°°(0, T; W), (2.3) 
p G H\<d,T;H) nC°([0,T]; V) HL 2 (0,T; W) H L°°(Q). (2.4) 

(ii) There are constants < p* < p* < 1 , p* > , and i^* > , depending 
only on the data, such that for every u E U a d the corresponding solution 
(p, p) satisfies 

< p* < p < p* < 1, < p < p*, a. e. in Q, (2.5) 

||p|||y 1 .°°(0,T;ir)nif 1 (0,T;V)nL oo (0,T;VK) 

+ ||/ i l| J H' 1 (o,r;H)nc ([o,T];y)nL 2 (o,T;vy)nL°°(Q) < (2-6) 

(hi) Let ux,v,2 E U a d, and let (pi,pi), (p 2 ,p 2 ) be the corresponding solutions 
to fll.2p - fll.5l) . Moreover, let u = ux — u 2 , p = pi — p 2 , P = Pi — P2 • TTien, 
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for all t G [0,T], 



max (||M*)|| 2 h + \\p(s)\\ 2 y) + f (||M*)H 2 v + M'Wh) ds 

0<s<t J Q 

< K*f\\u{s)f H ds, (2.7) 
Jo 

with a constant > that may depend on the data, but not on u\, U2 . 

Remarks: I. Owing to (12.71) . the solution operator S : u t-» (p, p) is 
Lipschitz continuous as a mapping from U a d (viewed as a subset of L 2 (Q)) 
into (^(0, T; H) n C°([0, T}; V)) x (L 2 (0, T; V) n C°([0, T]; #)) . 
2. Thanks to <^ and to / G C 2 (0, 1) , we have f'(p) G L°°(Q). More- 
over, owing to (12. 4p and to the embedding 1/ C L 6 (Q), we have p G 
C°([0, T]; L 6 (f2)) . Note that (2.3) implies, in particular, that p is contin- 
uous from [0, T] to H S (Q) for all s < 2. Now, provided that s is sufficiently 
large, we have H S (Q) C C(fi) ; consequently, p G C(Q) • Hence, without loss 
of generality (by possibly choosing a larger K* ), we may assume that also 

IHIc(Q) + IN|c°([0,T];L6(fi)) + ||p<||L 2 (0,T;L6(n)) < ^ ■ (2.8) 

We are now prepared to prove existence for the control problem (CP): 

Theorem 2.2 Suppose that the conditions (Al) and (A2) are satisfied. 
Then the problem (CP) has a solution u G U ad- 
Proof. Let {«„,} C {7 a d be a minimizing sequence for (CP), and let {(p n ,,p n )} 
be the sequence of the associated solutions to (ll.2IHjl.5p . We then can infer 
from (12. 6p the existence of a triple (u, p, p) such that, for a suitable subse- 
quence again indexed by n , we have 

u„ — y u weakly star in L°°(Q), 

Pn^p weakly star in W /1,oo (0, T; H) fl ^(0, T; V) D L°°(0, T; W), 
p n -> p weakly star in H l (0, T; H) n L°°([0, T]; V) l~l L 2 (0, T; W). 

Clearly, we have that u G U a d- Moreover, by virtue of the Aubin-Lions 
lemma (cf. [14, Thm. 5.1, p. 58]) and similar compactness results (cf. [16, 
Sect. 8, Cor. 4]), we also have the strong convergences 

p n ^p strongly in C°([0,T]; H s (n)) for all s < 2, (2.9) 
p n ^p strongly in C°([0,T]; H) D L 2 (0,T;V). (2.10) 
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From this we infer, possibly selecting another subsequence again indexed by 
n , that p n — > p pointwise a. e. in Q . In particular, p* < p < p* a. e. in Q 
and, since / G C 2 (0, 1) , also f'(p n ) — > f'(p) strongly in L 2 (Q) . Now notice 
that the above convergences imply, in particular, that 

Pn^p strongly in C°([0, T}; L 6 (fi), d t p n -> d t p weakly in L 2 (0, T; L 4 (fi)), 
p n — » p strongly in L 2 (0, T; L 4 (f2)), <9i/i n — >■ <9 t p weakly in L 2 (Q). 

From this, it is easily verified that 

Pn d t p n -> p <%p weakly in L\0, T; H), 
pndtPn P^p weakly in L 2 (0,T; L 3/2 (fi)). 

In summary, if we pass to the limit as n — > oo in the state equations fll.2[) — 
(11.51) written for the triple (u n , p n , p n ) , we find that (p, p) = S(u), that is, 
the triple (u, p, p) is admissible for the control problem (CP). From the weak 
sequential lower semicontinuity of the cost functional J it finally follows that 
u, together with (p, p) = S(u) , is a solution to (CP). This concludes the 
proof. B 

Remarks: 3. It can be shown that this existence result holds for much 
more general cost functionals. All we need is that J enjoy appropriate weak 
sequential lower semicontinuity properties that match the above weak con- 
vergences. 

4. Since the state component p is continuous on Q , the existence result 
remains valid if suitable pointwise state constraints for p are added (provided 
the admissible set is not empty). For instance, consider the case when the 
state has to obey the one-sided obstacle condition 

p{x\t)> 1 - VtG[0,T] (2.11) 

for some fixed x* G Vt (which of course requires that po{x*) > 1/2). If the 
set of admissible controls in \J a d is not empty, i. e., if there is at least one u G 
U a d such that the corresponding state component p satisfies (12. lip , then an 
optimal control exists. Indeed, we pick a minimizing sequence of admissible 
controls {u n } C U a d with associated states (p n ,Pn) obeying (12.111) . Since fl 
is a bounded three-dimensional domain, we can infer from (12. 9 p and Sobolev 
embeddings that p n — > p uniformly in Q , whence it follows that p satisfies 
fl2T[]). 
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3 Necessary optimality conditions 



In this section, we derive the first-order necessary conditions of optimality 
for problem (CP). In this whole section, we generally assume that the hy- 
potheses (Al) and (A2) are satisfied and that u £ U a d is an optimal control 
with associated state (p, p), which has the properties (I2.3l) - (l2.6p and (12.81) ; 
in particular, we have p £ C (Q) . For technical reasons, we need to take a 
slightly smoother nonlinear term / ; precisely, we take 

(A3) /£C 3 (0,1). 

With this assumption, we can improve the stability estimate (|2.7j) . Before 
stating the result, let us observe that (I2.3P implies, in particular, that the so- 
lution component p is weakly continuous from [0, T] into W , which justifies 
the formulation of the next estimate (13. ip . 



Lemma 3.1 Suppose that (A1)-(A3) are satisfied, and let Ui,u 2 £ U a d 
be given and (p 1; pi), (p 2 , P2) be the corresponding solutions to (jl.2p - fll.5D . 
Moreover, let u = Ui—u 2 , p = pi — p 2 , p = p\ — p% . Then, for all t £ [0, T] , 

max (\\ Pt ( S )\\l + \\p(s)r v + \\p(s)\\ 2 w ) + j\\\p t {s)\\ 2 H + \\p t (s)\\ 2 w ) ds 

< K* [ \\u(s)\\ 2 H ds, (3.1) 
Jo 

with a constant K% > that may depend on the data, but not on u\ , U2 ■ 
Proof. Obviously, the pair (p, p) is a solution to the system 

(e + 2p 1 )p t + 2pp 2 ,< + PPi,t + P2Pt- Ap = u a.e. in Q, (3.2) 

5p t -Ap = p- (/'(pi) - f{p2)) a. e. in Q, (3.3) 
dp dp 

— — = — = a. e. on S, (3.4) 
on oti 

p(x, 0) = p(x, 0) = , a.e. in Q. (3.5) 

We test Eq. ( 13. 2 p by p t . It then follows, with the use of Young's inequality, 
that 



[ WMWlds + i||VAi(*)HSr < 11 Hs)\\ 2 H ds 

(2\p\\p 2 ,t\ + \p\ \pi,t\ + \p2\\pt\)\pt\dxds . (3.6) 



./n 
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We estimate the terms on the right-hand side individually. In this process, 
Cj (i G IN) denote positive constants that only depend on the constants 
s, 5, p*, p*, p*, T, Kl, K% . On using Holder's and Young's inequalities, as well 
as the continuity of the embedding H 1 ^) C L 4 (f2) , we have that 



p| |pi, t | \nt\ dx ds < / \\n t (s)\\H\\n(s)\\ L <t(n)\\pi,t( s )\\LHn)ds 
Jo 

<7 fWmM da + ^ [\\piAs)\\vMs)\\vds. (3.7) 



o </n 



7 Jo 

Observe that, owing to (12.61) . the function s \-> \\pi,t(s)\\y belongs to L 1 (0, T) 
Moreover, by (I2.5P and Young's inequality, 

pt Q l>t 

\Pt\ \lh\ dxds < j \\Pt(s)\\ 2 H ds H \\p t (s)\ 2 H ds , (3.8) 



o Jn Jo 7 Jo 

where the second integral on the right-hand side can be estimated using (12. 7p . 
In addition, we have 

h := / 2 |p| \jjL 2 ,t\ \lh\ dxds 



'o Jn 



Now, observe that the embedding W C L°° (f2) , in combination with stan- 
dard elliptic estimates and (12.71) . implies that 

IIpOOIIl-M < ^ f \Hr)\\ 2 H dr + C 5 ||Ap(s)||i, (3.10) 
Jo 



whence 



h < 7 / IhOOI&ds 
Jo 

/lluwui,* + /W«MIIh ||a^)£a) , (3.11; 



7 Vo 

where the mapping s H- ||p2,i(s)||# belongs to L 1 (0,T). 

Next, we formally test Eq. (13. 3p by —Ap t . On integrating by parts, we find 

that 

$ f\Wpt{s)\\ 2 H ds + \\\kp{t)\\l < [ [(- fl + (f'( Pl )-f(p 2 )) A Pt dxds. 
Jo z Jo Jn 

(3.12) 
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After a further integration by parts, this time with respect to t , and invoking 
Young's inequality and (12. 7p . we find that 



t r r rt 



pAp t dxds = - I p(t) Ap(t) dx + I I p t Apdx ds < - \\ Ap(t)\\ 2 H 



>oJn Jn JoJn 

-t n rt 



+ 7 / \\fi t (s)\\ 2 H ds + ^ [ \\Ap(s)\\ 2 H ds + C 8 [ \\u(s)\\ 2 H ds. (3.13) 
Jo 7 Jo Jo 

Moreover, integration by parts with respect to t , with the help of Young's 
inequality, (12. 6p . and (12 .7p . yields: 

r (f'(p l )-f(p2))Ap t dxds < [ \f'(p 1 (t))-f(p 2 (t))\\Ap(t)\dx 



'oJn Jn 

+ ff (ir(Pl)l \Pt\ + \f"(pi) - f"( P 2)\ \ P 2,t\) \^ P \dxds 

Jo Jn 

< l\\Mt)\\ 2 H + j)\^p(s)\\ 2 H ds + C g j)\u(s)\\ 2 H ds + J 2 , (3.14) 
where 

h:= f (\f"{p x )- f"{p 2 )\\p 2 , t \\Ap\dxds. 



Jo Jn 

From this inequality, on applying the mean value theorem, (12. 6p . (12.71) . 
Young's inequality, and the continuity of the embedding if 1 (f2) C L 4 (f2) , 
we deduce the following estimate: 

h < C 10 [ \\Ap(s)\\ H \\p(s)\\ LHn) \\p2,t(s)\\ L 4(a) ds 
Jo 

< Cu max ||p(s)|| y / ||p2,t(s)||y \\Ap(s)\\ H ds 

< C 12 (J)\p 2 , t (s)\\ 2 v \\^p(s)\\ 2 H ds + j)\u(s)\\ 2 H ds) . (3.15) 

Observe that by (12. 6p the mapping s \-> \\p 2 ,t( s )\\v belongs to L l {Q,T). 
At this point, we may combine the estimates (I3.6p - ()3.15p : in fact, choosing 
7 > appropriately small, and invoking Gronwall's lemma, we find the 
estimate: 

*'lMs)lll + \\V p t {s)f H ) ds + max (||^( S )||^ + ||p( a )||^) 

/ 0<s<t 

< Cia f\Hs)\\ 2 H ds, (3.16) 
Jo 
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for all t G [0,T]. Next, we formally differentiate Eq. ( 13. 3 p with respect to 
t , and obtain 

5p tt - Apt = /it - /"(pi) Pt - (f"( Pl ) - f" ( P2 )) p 2>t , (3.17) 

with zero initial and Neumann boundary conditions for p t (cf. (13.4)) . (13.51) . 
and (11.51) ). Hence, testing (13.171) by p t , invoking Young's inequality, and 
recalling (j2.7[) . and (13 . 16)) . we find that 

+ j]\Vp t (s)\&ds < C lA j\\u{s)f H ds 

+ // \ P 2,t\\f"{pi)-f"{p2)\\pt\dxds. (3.18) 
Jo Jn 

Moreover, using Holder's and Young's inequalities, (A3), (12. 6p . and (12. 7p . we 
see that 



\p2 ! t\\f"(Pl)-f"(P2)\\Pt\dxds 

'o Jn 

ft 



< Ci 5 / ||p 2 ,t(s)|U4 (n) ||p(s)|| L 4 (n) ||pi(s)|| H ds 







2 J„ , ™„„ II ^„M|2 / || „ { „\ 1 1 2 



< C 16 / UptWII^da + max||p( S )||^ / U/MW^ 



o Jo 

t 



< C l7 / \\u{s)\\ l H ds. (3.19) 



Finally, we test (13. 1 7j) by —Ap t . Using Young's inequality and (I3.16p . we 
find that 



\\\V Pt {t)\\ 2 H + J \\^Pt(s)\\ 2 H ds < 7 ^ \\Ap t (s)f H ds 

+ — [ t \Hs)\\ 2 H ds+ ff \p2, t \\f"(p 1 )-f"( P 2)\\Ap t \dxds 

7 Jo Jo Jn 

< 2 1 f \\Ap t {s)\\ 2 H ds 

Jo 

+ — ( f\Hs)\\%da + max llpOOIfc fl^Wv ds) 
7 v Jo Jo y 

< 2 7 fllAp^Wlds + 9° f\Hs)\\ 2 H ds. (3.20) 

Jo 7 Jo 
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Choosing 7 > appropriately small, we can infer that the estimate (13. ip is 
in fact true. This concludes the proof. B 

3.1 The linearized system 

Suppose that h G L°°(Q) is an admissible variation with respect to u, i.e., 
that there exists A > such that u + Xh G U a d whenever < A < A . 
We have to determine the directional derivative DJ ve &(u)h of the "reduced" 
cost functional J re d{u) '■— J(u,Su) at u in the direction h. This requires 
to find the directional derivative DS(u)h of the solution operator S at u 
in the direction h. To this end, we consider the following system, which is 
obtained by linearizing the system (ll.2l) - (ll.5p at (p, p) : 

(E + 2p)r) t -Ar) + 2ji t £ + fl£ t + p t ri = h a.e. in Q, (3.21) 
£&- A£ = -f{p)i + 77 a.e. in Q, (3.22) 

^ = = a.e. on E, (3.23) 
on on 

£(x, 0) = rj(x, 0) = a.e. in a (3.24) 

We expect that (£,77) = DS(u)h, provided that fl3~211 - fl3~24T ) admits a 
unique solution (£,77). In view of (12.31) and (12.41) . we can guess the regu- 
larity of £ and 77 : 

£ G H\0, T; H) n C°([0, T]; V) n L 2 (0, T; W) PI L°°(Q), (3.25) 
77 G # x (0, T; if) n C°([0, T]; V) D L 2 (0, T; W). (3.26) 

Indeed, if (13.251) and (13.261) hold then the collection of source terms in (13.211) . 
i.e., the part h — 2p 4 £ — p£t — PtV , belongs to L 2 {Q) (as it should for a 
solution 77 satisfying (I3.26P ). whereas the regularity (13.261) for 77 allows us 
to conclude from (I3.22p that also £ G C(Q) (by applying maximal parabolic 
regularity theory, see, e.g., [7, Thm. 6.8] or [17, Lemma 7.12]). 

In fact, as to £, we can count on an even better regularity. Indeed, we may 
differentiate (I3.22p with respect to t to find that 

6^ - A£ t = -f (p) p t i - f'{p) & + m , (3.27) 

with zero initial and Neumann boundary conditions for £ 4 . Since the right- 
hand side of (I3.27P belongs to L 2 (Q) , we may test by any of the functions 
£t , £tt , and — A£ t , to obtain that even 

£ G # 2 (0, T; H) n C^flO, T]; y) n ^(0, T; W) . (3.28) 
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Notice, however, that this fact has no bearing on the regularity of 77, since 
the coefficient p t in (13.211) only belongs to L 2 (Q) . 

We first prove the well-posedness of the linear system fl3.2ip - fl3.24p . 

Proposition 3.2 Suppose that (Al)-(A3) are fulfilled. Then the system 
fl3.2ip - fl3.24p has a unique solution (£,77) satisfying (13.261) and f)3.28p . 

Proof. We proceed in series of steps. 

Step 1: Approximation. Following the lines of our approach in [3J, we use 
an approximation technique based on a delay in the right-hand side of f!3.22p . 
To this end, we define for r > the translation operator T T : L 1 (0, T;if) 
—> L 1 (0,T;if) by putting, for every v G L 1 (0,T;if) and almost every t G 
(0,T), 

(T T v)(t) =v{t-r) if t >t, and (%v) (t) = if t < r. (3.29) 

Notice that, for any t> G L 2 (Q) and any r > , we obviously have ||7^||l 2 (Q) 
< IMU 2 (Q)- 

Then, for any fixed r > , we look for functions (£ r , 77 r ) , which satisfy ( I3.25j) 
and (13.261) and the system: 

(e + 2p) rf t - Arj T + 2 fi t C + /j £[ + p^ T = /i a. e. in Q, (3.30) 
^[ - AC + f"(p) C = %rf a. e. in Q, (3.31) 

^- = ^- = a.e. on S, (3.32) 
an an 

C(x, 0) = r] T (x, 0) = a. e. in ft. (3.33) 

Precisely, we choose for r > the discrete values r = T/N , where A G IN is 
arbitrary, and put t n = nr , 0<?i<A, and J n = (0,t n ) . For 1 < n < N , 
we solve the problem 

(e + 2p) r) n j - At] n + 2fi t £ n + fi + p t rj n = h a. e. in ft x /„, (3.34) 
= a.e. on T x I n , r) n (x,0) = a.e. in ft, (3.35) 

an 

5 Cn,t ~ A£„, + f"(p) £ n = T T r] n a. e. in ft x /„, (3.36) 

din 

— — = a. e. on T x J n , £n(^, 0) = a. e. in ft, (3.37) 
an 

where the variables r] n and £ n , defined on /„ , have obvious meaning. Here, 
T T acts on functions that are not defined on the entire interval (0,T) ; how- 
ever, for n > 1 it is still defined by (13.291) . while for n — 1 we simply put 
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T T rj n = 0. Notice that whenever the pairs (£k,Vk) with 



£ fc G H\l k ; H) n C°(/ fe ; V) PI L 2 (/ fc ; fl C(Q x 4), (3.38) 
T] k G ^(4; tf) n C°(/ fe ; V) n L 2 (4; W), (3.39) 

have been constructed for 1 < k < n < N, then we look for the pair 
(£ n+ i,?7 n+ i) that coincides with ((, n ,Vn) in In, and note that the linear 
parabolic problem (13.36 p . (I3.37P has a unique solution £ n+1 on Q x I n+1 
that satisfies (13.381) for k = n + 1. Inserting £ n+1 in (13.341) (where n is re- 
placed by n+1 ), we then find that the linear parabolic problem (I3.34p . (I3.35P 
admits a unique solution r] n+ i that fulfills (13.391) for k = n + 1. Hence, we 
conclude that (CO = (£n,Vn) satisfies fl3T3ni) - fl3T33D . and (I3T25D . fl3T26|) . 

Step 2: A priori estimates. We now prove a series of a priori estimates 
for the functions (£ T ,?7 r ). In the following, we denote by Cj (2 G IN) some 
generic positive constants, which may depend on e, 5, p*, p*, p*, T, K^, , 
but not on r (i. e., not on N). For the sake of simplicity, we omit the super- 
script t and simply write (£, 77) . We recall the continuity of the embedding 
H^Q) C L 6 (tt). 

First a priori estimate. Observe that 2 p 77 774 = {pri 2 ) t — p t rf . Hence, 
testing (I3.30P by 77 , we have, for < t < T , 

J q (f + P) V(t) 2 dx + £ HV^II^da < i jT {\\v(s)\\ 2 H + \\h(s)\\%) ds 

+ 2 / / \pZ t \ |£| |t7| ds + / / \£ t \ \rj\ dx ds . (3.40) 
Jo Jn Jo Jn 

For any 7 > , we have by Young's inequality that 

\p,\\£ t \\r]\dxds < \\p,\\ L °°(Q) / \\ri(s)\\ H \\£ t (s)\\ H ds 

Jo 

<lj\Us)\\ 2 H ds + ^ fluids. (3.41) 

Moreover, 



Jn 



7 Jo 



Jn 



Ht\ |f| \n\dxds < / \\p t (s)\\ H ||f (s)|| L 4 (n ) lh(s)IU 4 (n) ds 



2 j„ , C'a 



< 7 / hOOII^* + — / ll/^)M£(*)||^- (3-42) 

Jo 7 Jo 

Notice that, by virtue of (12.61) . the mapping s !->■ belongs to 
LH0,T). 
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Next, we add £ on both sides of Eq. (I3.3ip and test the resulting equation 
by . On using Young's inequality again, we obtain: 

{ f \Us)f H ds + \(\\m\\ 2 H + \mt)\\ 2 H ) 



t r t 



< C 3 (y o \\v(s)\\ 2 H ds + l U(s)\\lds). (3.43) 

Adding the inequalities (I3.40p and (13.431) . and choosing 7 > sufficiently 
small, we conclude from the above estimates and Gronwall's lemma that 



(\ut)\\ 2 H + ut)\\v)dt + ™*(\m)\\v + mm 

<C 4 f \\h{t)f H dt. (3.44) 
Jo 

By comparison in (13.311) . and thanks to (13.321) . we may also infer (possibly 
by choosing a larger C 4 ) that 

T U(t)\\wdt < C A f \\h{t)f H dt. (3.45) 
Jo 



Second a priori estimate. We test (I3.30P by r\ t and apply Young's inequal- 
ity in order to obtain 

y\h(s)\\ 2 H ds + ~ II V V (t) \\ 2 H 

< 7T / \\ h (s)\\ 2 H ds+ [ [ (2\jk\\Z\ + \fi\\&\ + \to\\v\) \m\dxds. (3.46) 
Z£ Jo Jo Jn 

Since ft G L°°(Q) , we can infer from Young's inequality that 

161 \Vt\dxds < 7 / H(s)\\ 2 H ds + ^ [ UMfkds. (3.47) 



'0 Jn Jo 7 Jo 

Moreover, by virtue of Holder's and Young's inequalities, 

t r rt n rt 

pt\ \v\ \vt 



Jn 

t 



<1 



dxds<-f [ \\Vt(s)\\ 2 H ds + — [ \\pt(s)\\ 2 L 4 {n) \\r](s)\\ 2 LHn) ds 
Jo 7 Jo 

fht{s)\\ 2 H ds + °1 f \\pt{s)\\lUs)\\lds. (3.48) 
Jo 7 Jo 



17 



Observe that by (12. 6p the mapping s \-> \\pt(s)\\ v belongs to L 1 (0,T). 
Finally, we have, owing to the continuity of the embedding W C L°°(Q) and 

dsaiD, 

2\jh\ \i\W\dxds 




o Jn 



t m rt 



<lJ Q H(s)\\ 2 H d8+ ^ J \Ms)\\ 2 H U(s)\\U (n) ds 

< lfht{s)f H ds + ^{j\h{8)f B da + j]Ms)\\l\ms)\\ 2 H ds), 

(3.49) 

where, owing to (I2.6p . the mapping s (-)■ ||/ii(s)||# belongs to L 1 (0,T). 
Next, we test (13.311) by — A£ t to obtain, for every t e [0, T] , 

sfm t (8)\\%da + lm(t)\\ a B = ff (-(T T r)) + f"(p)0 A^dxds. (3.50) 
Jo / Jo Jn 

Now, by virtue of (I3.44p and (13.451) . and invoking Young's inequality, we have 

n(T T r]) A^ t dx ds 

< [\(%v)(t)\ \A£(t)\dx+ f I \d t {T T r))\ \Ai\dxds 

(3.51) 



o Jn 



< I WmWn + 1 j\r}t{s)f H ds + C 10 (l + I) f\h{s)f H ds. 



Moreover, it turns out that 



' r\p)iAi t dxds < i\f"(p(t))\\m\\m)\dx 

o Jn Jn 



+ 



f f \f"'(p)ptC + f"(p)Ct\ mdxds. (3.52) 
Jo Jn 

We have, owing to (12~5D and (13T44D . 

£|A^))||£(i)||A£(i)|d* < l\\Am\H + Cii /VoOII***- (3-53) 
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Also the second integral on the right-hand side of (13.521) is bounded, since 
fl?3j) . ( jggp , (EH, and dS35D imply that 

If \r\p)pti+r\p)it\ mdxds 

Jo Jn 

< C 12 / \\pt(s)\\ u(s)\\i m da+ / \ms)\\ 2 H ds 

Jo Jo 

< C 13 max ||e(*)||^ f)\pt(s)\\ 2 v ds + f]\Aas)\\ 2 H ds 

< C 14 f T \\h(s)\\ 2 H ds, (3.54) 



thanks to the continuity of the embedding 1/ C L 4 (f2). Thus, combining 
the estimates (I3.46P — f!3.54j) . choosing 7 > sufficiently small, and invoking 
Gronwall's inequality, we can infer that 

\\\Vt(t)\\ 2 H + \Ut)\\ 2 v)dt + max (\\ v (t)\\* + 

< C 15 f \\h{t)\\ 2 H dt. (3.55) 

Now, testing (13.301) by — A77, and arguing as for fl3.46p - fl3.49p . we find that 

\Ht)\\wdt < Cie f \\h{t)f H dt. 
Jo 

Next, we differentiate Eq. (13.311) with respect to t. We obtain: 

8t tt -A{, t = d t {T T r))-f"{p)ptt-f"{p)tit a.e. in Q. (3.56) 

From (I2.5p . (12.61) . (13.441) . and (I3.55p . we can infer that the expression on the 
right-hand side of (13.561) is bounded in L 2 {Q) . Therefore, we may test (I3.56P 
by any of the functions £ t , — A£ t , and C, u , m order to find that 

l T (Uu(t)\\ 2 H + m(t)\\ 2 H )dt + max ||6(t)|| 2 v, < C l7 £ \\h(t)\\ 2 H dt . 

(3.57) 

Step 3: Passage to the limit. Let (£,n,Vn) denote the solution to the sys- 
tem f l3.3()p -f l3T33|) associated with t n = T/N, for N G IN. In Step 2, we 
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have shown that there is some C > 0, which does not depend on N, such 
that 

Mn\\h* (0,T;/i")nC 1 ([0,T] ]V)nm (0,T;W)nC(Q) 
+ \\VN\\H 1 (0,T;H)r\C°{[0,T\;V)nL 2 (0,T;W) < C- (3.58) 

Hence, there is a subsequence, which is again indexed by N , such that 

6v-^£ weakly star in H 2 (0, T; H) fl W 1,oo (0, T; V) fl //' 1 (0, T; W), 
77^^77 weakly star in i7 1 (0, T; if) fl L°°(0, T; V) fl L 2 (0, T; W). (3.59) 

By compact embedding, we also have, in particular, 

£tv — > £ strongly in C(Q), t]n — > V strongly in L 2 {Q), (3.60) 

so that p rjN tt ~> PVt an d p6v,t ~^ p£t > both weakly in L 2 {Q) , /"(p) ^ — > 
/"(p) £ strongly in L 2 (Q) , as well as fit 6v,t -> fit & and fi t T] N ->> p t 77 , both 
strongly in L 1 (Q). Finally, it is easily verified that {Tt/nVn} converges 
strongly in L 2 (Q) to 77. In conclusion, we may pass to the limit as N — > 00 
in the system fl3.30l) - fl3.33l) (written for r = T/N) to find that the pair (£, 77) 
is in fact a strong solution to the linearized system fl3.2ip - fl3.24p . 

It remains to show the uniqueness. But if (£1, 7/1) , (£ 2 , ^2) are two solutions 
having the above properties, then the pair (£, 77) , where £ = £1 — £2 and 
77 = 771 — 772 , satisfies fl3.2ip - fl3.24p with h = . We thus may repeat the first 
a priori estimate in Step 2 to conclude that £ = 77 = . This concludes the 
proof. a 

3.2 Directional differentiability of the control-to-state 
mapping 

In this section, we prove the following result. 

Proposition 3.3 Suppose that the assumptions (Al)-(A3) are satisfied. 
Then the solution operator S , viewed as a mapping from U a d, subset of 
L 2 (Q) , into 

(H\0, T; H) n C°([0, T];V)n L 2 (0, T; Wj) x (C°([0, T]; if) D L 2 (0, T; V)) , 

is directionally differentiate at u in the direction h . The directional deriva- 
tive (£,77) = DS(u)h is given by the unique solution (£,77) to the linearized 
system Iff21j) - (I323)) . 
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Proof. Let A > be such that u + Xh G U ad for < A < A . We put 

u x = u + Xh, (p\/i A ) = S(u x ), y x = p x -p-\^ z x = p, x — p — Xrj. 

We have to show that there is a function Z : [0, A] — > [0, +oo) with 
liniA\o Z(X)/X 2 = such that 

lb A |l^ 1 (0,T;_H")nC ([0,T];V)nL 2 (0,T;W/) + II Z> ~ II C°([0,T] ;H )(1L 2 (0,T;V) — % (X) ■ (3.61) 

Using the state system (ll.2l) - (ll.5p and the linearized system (I3.21I) - (I3.24I) . 
we easily verify that for < A < A the pair (y x , z x ) is a strong solution to 
the system 

(e + 2p) z x + p t z x + fi y x + 2fi t y x - Az x 

= -2 (/i A - fi t ) (p x -p)- (p x - p t ) (p x - ft) a. e. in Q, (3.62) 
5y x - Ay x + /'(/) - f(p) - X f"(p) £ = z\ a. e. in Q, (3.63) 

— — = — — = (J, a. e. on L, (3.64) 
on an 

y x (x, 0) = z x (x, 0) = a. e. in fi. (3.65) 

Notice that 

y x G H 1 ^, T; H) PI C°([0, T]; V) D L 2 (0, T; W) PI C(Q), 
^ A G H\0, T; H) n C°([0, T];V)n L 2 (0, T; W). 

For the sake of a better readability, in the following estimates we omit the 
superscript A of y x and z x . As before, we denote by Ci (i G IN) certain 
positive constants that only depend on e, 5, p*, p*, /!*, T, X*, i^, -K3 , but not 
on A . 

We now add y on both sides of Eq. (I3.63P and test the resulting equation 
by y t . Using Young's inequality, we find that for all t G [0,T] it holds 

|jT||»t(«)llird»+i(l|Vy(t)|^+||y(t)||y < - 5 j\\z{s)f H ds 

+ Ci [ t \\y(s)\\ 2 H ds + C 2 [ t \\(f'(p x )-f'(p)-Xf"(p)0(s)\\ 2 H ds.(3.m) 
Jo Jo 

In order to handle the third term on the right-hand side of (I3.66p . we note 
that the stability estimate (13.11) implies, in particular, that 

||p A -p|||oc (Q) < K;X 2 \\h\\ 2 LHQ) , (3.67) 
21 



that is, p x — > p uniformly on Q as A \ 0. Since / G C 3 (0, 1) , we can infer 
from Taylor's theorem that 

\f'(p X ) ~ f'(p) ~ A/"(p)£| < \ max |/»| |p A - p\ 2 +\f"(p)\ \y\ on Q. 

(3.68) 

It then follows from the estimates (13.1 1) and (I3.66P that 

* + Jlll/(*)llv < 7 /ikOOII**' + C,j\y{s)f H ds + C 4 A 4 , 

(3.69) 



2 Jo 2 5 j 



since we can assume that ||/i||z,2(q) < 1- Next, observe that 2pzz t = 
(pz 2 ) t — p t z 2 . Therefore, testing (13.621) by z yields for every t G [0,T] that 

/ Or + Pt*)) 2 ^*)^ + / II V2 ( s )IIh^ = - / / (fiyt + 2p t y)zdxds 

-2 / / (/i A - /it) (p A - p) zdxds- / (p A - p t ) (/i A - /i) z cZs. 
io in io in 



(3.70) 

We estimate the terms on the right-hand side of (I3.70P individually. At first, 
using (12. 6p and Young's inequality, we find that 

t f pt pt 

\fi\\yt\\z\dxds < 7 / \\y t (s)\\ 2 H ds + / (3.71) 
o in io 7 io 

Moreover, using the continuity of the embedding if 1 (fi) C L 4 (Q) , as well as 
Holder's and Young's inequalities, 

t r rt 

\fk\ \y\ \A dxds < j \\ji t {s)\\ H \\z{s)\\ LHsi) \\y(s)\\ L A {n) ds 
'o Jn Jo 

ft /~t pt 

|2 - °< 



< 7 / ||*( S )||£<k + / ll^Wll^llyWII^a. (3.72) 
io 7 Jo 

Observe that by (12. 6p the mapping s >->■ ||/it(s)||jj belongs to L 1 (0,T). 
At this point, we can conclude from (13. ip and (I3.67p . invoking Young's in- 
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equality, that 

2 \jjL X — p t \ \p x — p\ \z\dxds 



'0 Jn 



< 2 / \\(K-fit)(s)\\ H \\(p-p)(s)\\ Loo{Q) \\z(s)\\ H ds 







< C 7 \\(P X ~ p)(s)\\ 2 Lao{Q) - fi t )(s)\\ 2 H ds + J \\z{s)f H ds 

< [ \\z(s)\\ 2 H ds + C 8 A 4 . (3.73) 
Jo 

Finally, we invoke (13. ip and Holder's and Young's inequalities, as well as the 
continuity of the embedding H 1 ^) C L 4 (f2) , to obtain that 

| p x — p t | | /i A — p I I z | dx ds 

o Jn 

< max ||z(s)||iy / \\(Pt ~ Pt)( s )\\ L 4 (n) \\(V X -fi)( s )\\ L i { n) ds 

— — J 

< 7 max ||z(s)||^ + — / \{p x t - p t )(s)\\ v ds / || (p x - p)(s) \\ v ds 



0<s<t 



7 Jo 



<7max||z( s )||^ + C 10 A 4 . (3.74) 

0<s<t 

Combining the estimates fl3.69p - fl3.74p . taking the maximum with respect 
to t G [0, T] , adjusting 7 > appropriately small, and invoking Gronwall's 
lemma, we arrive at the conclusion that (y x , z x ) = (y, z) satisfies the in- 
equality 

ll?/ A |lf/ 1 (0,T'; J H')nC ([0,T];y) + II Z II C°{[0,T] ;H)nL 2 {[0,T] ;V) ^ Cn \ A . (3.75) 

Finally, testing f)3.63p by —Ay x , and using f)3.68p . we find that also 

\\y X \\h(o,T;W) < C12A 4 . (3.76) 
This concludes the proof of the assertion. B 



Corollary 3.4 Let the assumptions (Al)-(A3) be fulfilled, and let u G U a d 
be an optimal control for the problem (CP) with associated state (p,p) = 
S(u) . Then, for every v G U a d , 

n(32u(v—u) dx dt+ (p(T)—p T )^(T)dx + / / Pi (ft— pr) V dx dt > 0, 
s Jn Jo Jn 

(3.77) 
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where (£, rf) is the unique solution to the linearized system (I3.21I) - (I3.24I) as- 
sociated with h = v — u . 

Proof. Let v G U a d be arbitrary. Then fa = v — u is an admissible 
direction, since u + Xh G f/ a( ^ for < A < 1 . For any such A , we have 

^ < J(u + Xh, Sju + Xh)) - J(u, S(u)) 

X 

J(u + Xh, S(u + Xh)) - J{u, S{u + Xh)) J{u, S(u + Xh)) - J(u, S(u)) 



< 



X X 



It follows immediately from the definition of the cost functional J that 
the first summand on the right-hand side of this inequality converges to 
f f n f3 2 uhdxdt as A \ . For the second summand, we obtain from 
Proposition 3.3 that 

, J(u,S(u + Xh)) - J(u,S(u)) 

lim 

a\o A 

= (p(T) - p T ) £(T) dx + / I fix {p,- ht) rjdxdt , 
Jn Jo Jn 

whence the assertion follows. ■ 



3.3 The optimality system 

Let u G U a d be an optimal control for (CP) with associated state (p, p) = 
S(u) . Then, for every v G U a d , (I3.77P holds. We now aim to eliminate 
(£> v) by introducing the adjoint state variables. To this end, we consider 
the adjoint system: 



- (e + 2p) q t - p t q - Aq = p + ft (/i - p T ) a. e. in Q, (3.78) 

— — = a. e. in E, q(x,T) = a. e. in fi, (3.79) 
on 

-5p t - Ap + f"(p)p = pq t ~ JkQ in Q, (3-80) 

|^ = on S, 5 p(T) = p(T) - p T mQ, (3.81) 

which is a linear backward-in-time parabolic system for the adjoint state 
variables p and q. 
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It must be expected that the adjoint state variables (p, q) be less regular 
than the state variables (p,p). Indeed, we only have p(T) G L 2 (Q) , and 
thus (13.801) and (I3.8ip should be interpreted in the ususal weak sense. That 
is, we look for a vector-valued function p G tf^T; V*) H C°([0, T]; H) n 
L 2 (0, T; V) that, in addition to the final time condition (I3.8ip . satisfies 

(-6 Pt {t),v) V *y + [ Vp{t)-Vvdx+ [ f"{p{t))p{t)vdx 

Jn Jn 

= [ (p.(t)qt(t)-th(t)q(t)) vdx, (3.82) 

Jn 

for every v G V and almost every t G (0, T) . Notice that if q G i? 1 (0, T; H)(l 
C°([0,T];V) , then it is easily seen that pq t — ptq G L 3 ^ 2 (Q), so that the 
integral on the right-hand side of (13.821) makes sense. On the other hand, if p 
has the expected regularity then the solution to (I3.78p . (I3.79P should belong 
to H\0,T;H) nC°([0,T];V) C)L 2 (0,T;W). 



Lemma 3.5 Suppose that the system (I3.78p - (l3.8ip has a unique solu- 
tion (p,q) where p G H\0,T;V*) PI C°([0, T]; H) D L 2 (0,T;V) and q G 
H 1 (0,T;H)nC ([0,T};V)nL 2 (0,T;W). Then we have 

(p(x,T) — Pt(x)) £(x,T) dx + / / f3± (p, — /It) V dx dt — / / qhdxdt. 

Jo Jn Jo Jn 

(3.83) 

Proof. The assertion follows from repeated integration by parts, using the 
well-known integration by parts formula 

((V t (t),w(t)) v *y + (W t (t),v(t)) v *y)dt= [ (v(T)w(T) - v(0)w(0)) dx, 

Jn 

which holds for all functions v,w G iJ 1 (0,T;l / *) H L 2 (0,T;V) . Since this 
calculation is standard in optimal control theory, we may leave it to the 
reader to work out the details. ■ 



Proposition 3.6 The adjoint system (I3.78p -f l3.8ip has a unique solu- 
tion (p,q) with p G H l {Q,T-V*) n C°([0,T);H) n L 2 (0,T;V) and q G 
H\0, T; H) n C°([0, T];V)n L 2 (0, T; W) , where fl3T80|) and fl3T8T| l are un- 
derstood in the sense of ( I3.82[) . 

Proof. We proceed in a series of steps. 
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Step 1: Approximation. As in the proof of Proposition 3.2, we employ a 
delay technique. However, this time we have to use a negative delay, since 
the system (13.78p - (13.8ip runs backwards in time. 

For < t < T, we define the translation operator T T : L 1 (0,T;if) — > 
L 1 (0,T; H) by putting, for every v G L 1 (0,T; H) and almost all t G (0, T) , 



(f r v) (t) =v(t + r) if t < T — t, and (%v) (t) = if t > T — r; 

(3.84) 

clearly, we have that ||7t^||z, 2 (q) < IM|l 2 (Q), for every -u G L 2 (Q) and any 
r G (0,T). 

We then consider for r G (0, T) the following approximating problem: find 
functions 

p T G /^(O.T; V*) nC°([0,T];#) n L 2 (0, T; 1/), 

g r G ^(0, T; if) n C°([0, T]; V) fl L 2 (0, T; VP) (3.85) 

that solve the system 

-(e + 2p) q T t - Aq r + q T = (1 + p t ) f T q T + t T p T + Pi(fi - #r) 

a.e. in Q, (3.86) 

= a.e. on S, q T (x,T) = a.e. in f2, (3.87) 

an 

(-5 P T t (t),v) V *y+ [Vp r (t)-Vvdx + [f"(p(t)) P T (t)vdx 

Jn Jn 
= J(p{t)q T t {t) - p t {t)q T {t^v dx \/vEV, for a.e. t G (0,T), (3.88) 

5p r (T) = p(T)-pr a.e. in fi. (3.89) 

We choose for r G (0, T) the discrete values r = T/N, where A G IN 
is arbitrary, and we put t n — nr, < n < N , and /„, = (t n ,T) . For 
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n = N — 1,... ,1,0, we solve the problem: 

-(e + 2p) q n>t - Aq n + q n = (1 + p t ) f T q n + %p n + fii(p> - H>t) 

a.e. in Q x I n , (3.90) 

0q n 

-— ^ = a.e. on S, q n (x,T) = a.e. in f2, (3.91) 
on 

(-5p n ,t(t),v) V * >V + Vp n {t) ■ Vvdx + f\p(t))p n (t)vdx 

Jn Jn 
= J (jl(t) q n>t (t) - pt(t) q n (t) j u cfe for all v e V, for a. e. t G (t n , T), 

(3.92) 

5p n (T) = p(T) — pt a.e. in f2. (3.93) 

Here, T T acts on functions that are not defined on the entire interval (0, T) ; 
however, for n < N — 1 it is still defined by (I3.84p . while for n = N — 1 we 
simply put 7~ T p n = T T q n = 0. Whenever the pairs (pk,qk) with 

Pk e H\l k ; v*) n C°(4; if) n L 2 (I k ; V), (3.94) 
g fc G H\l k ; H) H C°(7 fc ; V) fl L 2 (J fc ; W), (3.95) 

have been constructed for < n < k < N — 1 , we look for the pair 
(£n-i> ^n-i) that coincides with (£ n , rj n ) in /„ . Note that the linear parabolic 
problem (I3.90p . (13.911) has a unique solution q n _i on Q x I n _ x that satisfies 
(I3.95P for k — n — 1 (see, e.g., [12]) . On inserting g n _i in (13.921) (with 
n replaced by n — 1), we then find (e.g., by using an appropriate Galerkin 
approximation) that the linear parabolic problem (13. 92p . (I3.93P admits a 
unique solution p n -\ that fulfills (I3.94p for k = n — 1 . Hence, we conclude 
that (p T ,q T ) = (p ,9o) satisfies fl3TS5|l - (l3Tgg|l and f[3TS5]) . 

Step 2: A priori estimates. We now prove a series of a priori estimates 
for the functions (p T ,q T ). In the following, we denote by C; (i G IN) some 
generic positive constants, which may depend on e, 5, p*, p*, /!*, T, K^, K\ , 
but not on r (i.e., not on N). For the sake of simplicity, we omit the 
superscript r and simply write (p, q) . 

We multiply (I3.86P by —q t and integrate over Q x [t, T] to obtain, using 
Young's inequality, 

\ [ T \Ms)\\ 2 H ds+ 1 -(\\Vq(t)\\H+ \\q(t)\\%) 



< 



C! + C 2 [ + Il?(*)ll2r)<te + / / li^l 1^1 \Qt\dxds. (3.96) 
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Moreover, by virtue of Holder's and Young's inequalities, and invoking the 
continuity of the embedding H 1 ^) C L 4 (f2) , we have, for any 7 > 0, that 

n\pt\\t T q\\q t \dxds < / \\p t (s)\\ L 4 {n) \\T T q(s)\\ L 4 {n) \\q t (s)\\ H ds 
A Jt 

< 7 f ht{s)\\ 2 H ds + ^ f T \\p t {s)\\l\\q{s + r)\\ 2 v ds. (3.97) 
Jt 7 Jt 

Observe that by (12. 6p the mapping s h-> ||pt(s)||y belongs to L 1 (0, T). 

Next, we insert v = p(t) in (13.881) and integrate over [t, T] , where t £ [0, T] . 
We find, using f!3T9"3D . that 

~lb(*)ll* + \[ \\Vp{s)\\ 2 H ds < C 4 + C 5 £ T \\p(s)\\ 2 H ds 

+ / / Ia*I IpI ^ x + / / I A** I l?l IpI dx ds . (3.98) 
7t Jn Jt Jn 

Let us denote for short the last two integrals on the right-hand side of (I3.98|) 
by Ii and J 2 , respectively. Since ft £ L°°(Q) , we have that 

Jx < 7 f T \\qt(s)\\ 2 H ds + ^ [ T \\p(s)r H ds. (3.99) 
7t 7 A 

Moreover, we conclude from (12.61) . invoking Holder's and Young's inequalities, 
that 

h < f \\ik{s)\\ H ||?(s)|U*(n) \\p{s)\W(si)ds 



t 

< 1 f\\p{s)\\ 2 v ds + ^ f \\Ms)\\ 2 H \\q(s)\\ 2 v ds, (3.100) 
Jt 7 Jt 

where the mapping s H- |bt( s )||i? belongs to L 1 (0,T). 
Now, we combine the estimates f )3.96j) (13. lUUp . On choosing 7 > suffi- 
ciently small, and on applying Gronwall's lemma, we find that 

\\p\\c°([o,T];H)nL 2 (o,T;V) + II q\\ H^T^ncoQa^v) < C 8 . (3.101) 

It is now a standard matter to verify, by comparison in (I3.86j) and (13.881) . 
respectively, that also 

IblU^O.TjV*) + |MU 2 (0,T;W) < Cg. (3.102) 
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Step 3: Passage to the limit. Let (p^, q^) denote the solution to the sys- 
tem (lii^JHIMZj) associated with t n = T/N, for N e IN. In Step 2, we 
have shown that there is some C > 0, which does not depend on N , such 
that 

lbiv||ff 1 (0,T;F*)nC°([0,T];J3")ni 2 (0,T;V) + \\QN l|// 1 (0,T;H)nC o ([0,T];V)nL 2 (0,T;iy) < C- 

Hence, there is a subsequence, which is again indexed by N , such that 
p N -> p, weakly star in H\0, T; V*) n L°°(0, T; if) n L 2 (0, T; 1/), 
g^v -> g, weakly star in if ^O, T; if) n L°°(0, T; 1/) D L 2 (0, T; W). 

By compact embedding, we also have the strong convergences (see, e. g., [IB]) 

Pn — ^ P, strongly in L 2 (Q), 
q N ^q, strongly in C°([0,T}; H) f] L 2 (0,T;V). 

From this, we can conclude the following convergences: 

PlN,t PQt weakly in L 2 (Q), Aq N — > Aq weakly in L 2 (Q), 
Tt/nQn -> q strongly in L 2 (Q), T t /nPn ->■ P strongly in L 2 (Q), 
Pt Tt/nQn -»■ PtQ strongly in L l (Q). 

Hence, passing to the limit as N — > oo in fl3.86p - fl3.87l) for r = T/N , we 
find that the pair fjo, q) gives a strong solution to the parabolic problem 
fl3.78p - fl3.79p . Next, we notice that the weak convergence of {pn} to p in 
H\Q,T\V*) C]L 2 (0,T;V) implies that p N -> p weakly in C°([0,T]; H) . 
We may thus conclude, in particular, that 5p(T) = p(T) — px ■ Since f"(p) 
and p are bounded, we also have the following convergences: 

f" (P)Pn -> f" (P) V strongly in L 2 (Q), p q N>t -> pq t weakly in L 2 (Q). 

Therefore, f)3.82p is fulfilled for any v in C 1 (Q) , which is a dense subset 
of H 1 (O) , because Q is a Lipschitz domain. From this, it easily follows 
that f)3.82p is satisfied for every v G ii 1 (fi). In conclusion, the pair (p,q) 
is a solution to the adjoint system f)3.78l) - fl3.81l) that enjoys the asserted 
smoothness properties. 

Uniqueness remains to be shown. But if (pi,qi), (p2> #2) are f wo solutions 
having the above properties, then the pair (p, q) , where p = p\ — p<i and 
q = q\ — (72, satisfies fl3.78p - fl3.8ip . where the inhomogeneities (3\{p — pr) 
and p — pr on the right-hand sides of (I3.78P and f)3.8ip , respectively, cancel 
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out by subtraction. We may then repeat the a priori estimates of Step 2 
to see that in the present situation the constants C\ and C4 appearing, 
respectively, in (I3.96P and (13.981) . simply do not occur. Consequently, the 
application of Gronwall's lemma yields p = q = . This concludes the proof. 

■ 

In summary, we have proved the following result concerning first-order nec- 
essary optimality conditions. 

Theorem 3.7 Suppose that u £ U a d is an optimal control for (CP) 
with associated state (p,fi) = S(u) . Then the adjoint system (I3.78j) - (I3.81|) 
has a unique weak solution (p,q) with p £ H l (§,T\V*) D C°([0,T]; H) H 
L 2 (0,T;^), q £ H 1 (0, T; H) H C°([0, T];V) D £ 2 (0, T; W 7 ) ; moreover, /or 
am/ v £ t/ad, we have the inequality: 



Remark: 5. Since U a d is a nonempty, closed and convex subset of L 2 (Q) , 
(13.1031) has the following implications: 

• For 02 > 0, the optimal control u is nothing but the L 2 (Q) orthogonal 
projection of — Z?^ 1 Q onto U a d ■ In other words, 

u(x,t) =V[o,u(x,t)] {-P2 1 lix^)) a - e - in Q, 

where, for any a, b £ 1R such that a < b, 



a bang-bang situation where no information can be recovered for points 
at which q(x,t) = 0. 
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(3.103) 



'P[a,b](^) := min {6, max{a,«}}. 



• For 02 = , we have that almost everywhere 

f if q(x,t) > 
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